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Abstract. A search for gamma-ray bursts (GRBs) above 
20 TeV within the field of view (1 sr) of the HEGRA 
AIROBICC Cherenkov array (29°N, 18° W, 2200 m a.s.l.) 
has been performed using data taken between March 1992 
and March 1993. The search is based on an all-sky survey 
using four time scales, 10 seconds, 1 minute, 4 minutes 
and 1 hour. No evidence for TeV-emission has been found 
for the data sample. Flux upper limits are given. A special 
analysis has been performed for GRBs detected by BATSE 
and WATCH. Two partially and two fully contained GRBs 
in our field of view (FOV) were studied. For GRB 920925c 
which was fully contained in our FOV, the most signifi- 
cant excess has a probability of 7.7-10 -8 (corresponding 
to 5.4cr) of being caused by a background fluctuation. Cor- 
recting this probability with the appropriate trial factor, 
yields a 99.7% confidence level (CL) for this excess to be 
related to the GRB (corresponding to 2.7<r). This result is 
discussed within the framework of the WATCH detection. 

Key words: gamma rays: bursts - gamma rays: observa- 
tions - cosmic rays 



1. Introduction 

Emission of TeV/PeV gamma-rays associated with GRBs 
has been extensively searched. These studies are moti- 
vated by both the experimental results obtained with the 
satellite experiments as e.g. EGRET and by theoretical 
models (Meszaros et al. 1994) which predict or at least 
allow TeV emission to be produced in GRBs. So far none 



of these searches have revealed any convincing evidence 
for VHE emission (see for example Aglictta ct al 



1993 



Alexandreas et al. 1994; Borione et al. 1995; Connaughton 
et al. 1997] ; Dazeley et al. 1997), although some tenta- 
tive positive evidence has been found recently in this en- 
ergy range (Plunkett et al. 1995] ; Krawczynski et al. 1995; 



Amenomori et al. 1996|) 



Up to now five BATSE GRBs have been detected by 
EGRET with p hotons of energies up to 18 GeV (Hur- 
ley et al. 1994[ ). Those GRBs are among the most in- 



tense ones recorded by BATSE in the FOV of EGRET, 
so observations are compatible with the hypothesis that 
all GRBs emit GeV photons but only the strongest ones 
are above the EGRET sensitivity. A simple power law ex- 
trapolation of the Superbowl GRB spectrum (Sommer et 
al. [1994D predicts that ~20 photons above 20 TeV should 
be observed with the extensive air shower (EAS) array of 
wide angle integrating Cherenkov counters (AIROBICC) 
within 25 seconds while expecting only 0.1 background 
events. This would lead to a highly significant detection by 
several other experiments currently operating (Cherenkov 
telescopes, EAS arrays). The previous extrapolation ne- 
glects source-intrinsic cutoffs and the attenuation through 
interaction with the low energy cosmic photon background 



(Wdowc zyk ct al. |1972| Mannheim et al. |1996| ; Stecker & 
de Jager 1997), which is expected for cosmological sources 
with rcdshift greater than 0.1. Another interesting feature 
of strong GRBs is that their high energy emission can be 
delayed and have l onger durations than keV-MeV emis- 
sion (Hurley et al. 1994 ). These ideas motivate searches 
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for counterparts at times and with durations indepen- 
dent of those given by the space detectors at lower en- 
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ergies, especially when considering that the latter could 
miss GRBs out of their FOV, or with a very hard spec- 
trum, as has been already suggested (Piran & Narayan 
1995; see also Kommers et al. 1997). Searches for GRBs 



with the HEGRA experiment using other data periods can 



be f ound in Krawczynski ( |1997D, Fu nk ( |1997| ), Pad illa et 
al. ( |1997| ), Krawczynski et al. J199SQ , Padilla ( |1998| ). 



2. The AIROBICC array 

The AIROBICC array is part of the cosmic ray (CR) 
HEGRA detector complex located at the Roque de los 
Muchachos on the Canary Island La Palma (28.75° N, 
17.89° W, 2240 m a.s.L). It is an array of 7x7 stations (cur- 
rently enlarged to 8x8 plus a subarray of 4x4 to increase 
detector density in the center) with a regular grid spacing 
of 30 m, thus covering more than 32000 m 2 . Each station 
consists of a 40 cm diameter reflecting cone which focuses 
the incoming light onto a fast 20 cm diameter photomulti- 
plier tube (PMT). The set is placed inside a protective hut 
with a lid which can be opened through remote control. 
The PMT is covered with a blue filter (A =300-480 nm) to 
improve signal to noise ratio (S/N). The PMT output is 
amplified in the hut and sent through a 150 m long cable 
to a constant fraction discriminator (CFD) placed in the 
electronics container. The CFD is set to a level equivalent 
to 5a of the night sky background fluctuations. Whenever 
six or more AIROBICC stations exceed the CFD thresh- 
old within 200 ns, a trigger is produced. Calibrations are 
performed every 20 minutes to measure the relative delay 
between stations, the response of the TDC and the ADC 
pedestals. The mean dead time after each recorded event 
is —8 ms. Detailed descriptions of AIROBICC and the 
other compon ents of the HEGRA c omplex can be found 
in Karle et al. 1995a, Fonseca et al. 1995a, Fonseca et al. 



1995b 



The ADC signals (AIROBICC + scintillator array) are 
used to locate the EAS core position with a typical error 
<18 m. The Cherenkov light front is fitted by a cone to 
determine the arrival direction of the incident CR. The 
angular resolution (angular distance containing 63% of 
events for a point source) is 0.29° when 12 stations are 
fired (standard cut) and it continuously improves with in- 
creasing number of triggered huts. The absolute pointing 
accuracy has been estimated through comparison with the 
first HEGRA Cherenkov telescope to be better than 0.2° 



(Karle et al. 1995a). The FOV, which is limited by the ac- 
ceptance of the Winston cone, is 1 sr and the trigger rate 
is —20 Hz. The energy threshold (50% trigger probability) 
is estimated to be -16(25) TeV for 7 s and -29(37) TeV 
for hadrons with zenith angle 6 < 20° (20° < 9 < 3 5°) us - 
ing flux and Monte Carlo studies (Martinez et al. 1995[ ). 
The price that has to be paid for the advantages of the 
Cherenkov technique (lower energy threshold, better an- 
gular resolution) is that observations are restricted to 
moonless clear nights (—10% duty cycle). 
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Fig. 1. Systematic error in the estimate of the background 
with the armulus method as a function of the ratio of the 
solid angles of the OFF and ON windows. The ON window 
is a fixed circle of 0.7° radius and the OFF window is 
an annulus of fixed 1.5° inner radius and variable outer 
radius. The time interval is 4 minutes which is expected to 
have the lowest statistical error. The points are the average 
over several hundred windows. The sample used for this 
calculation consists of MC events generated as described 
in the text. The systematic error has a constant value of 
less than 5% up to a certain size of the OFF window, from 
which it grows very fast. The arrow points to the value 
used in this work. This error has to be balanced against 
the statistical error which in the best case is 20%. Note 
that the sign of the systematic error is positive, i.e. it is 
conservative. 

3. Analysis 

The data analyzed in this work were taken during the first 
year of operation of AIROBICC, between March 1992 and 
March 1993. The mean trigger rate in this period was —16 
Hz. The sample contains — 4-10 7 events (—800 hours) and 
after a successful arrival direction reconstruction of those 
showers with >7 fired stations, about — 2.5 • 10 7 events 
(—60%) survive. It has been shown that due to the low 
counting statistics encountered when dealing with tran- 
sient phenomena, one has to keep as many events as pos- 
sible; and so we do not apply the standard cut in the num- 
ber of fired huts. The angular resolution is then — 0.6°. 
The analysis is based on a binned (in time and in space) 
all-sky search. We define an ON window or source window 
and an OFF window or background window where events 
are counted for a certain time interval. The search strategy 
considers each event as being in the middle (in time) of a 
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Fig. 2. Cumulative number of trials (windows) versus probability in the all-sky search. The time scale is indicated in 
brackets. The solid lines represent real data while the dashed lines represent the results (normalized to the real data) 
over a MC sample 10 (20 in the 1 hour scale) times larger than the real one. Deviations are not significant (p > 0.1). 



burst, so for every single event in the data sample we take 
a pair of ON-OFF windows with a time interval centered 
at the time of the event. The ON window is a circle of 0.7° 
radius centered at the position of the event and the OFF 
window is an annulus of 1.5° and 10.0° radii concentric 
to the ON window. The solid angles covered by the ON 
and OFF windows are denoted by Qon and Qoff respec- 
tively, and their ratio £Ion /^off by a. There is always 
one event inside the ON window and in order to not over- 
estimate the significance of a possible signal this event is 
not counted. The size of the ON window is chosen to max- 
imize the S/N according to our resolu tion (see above) as 
described in Alexandreas et al. ( |l993a ) and it is expected 
to contain 75% of the events for a point source. While the 
size of the OFF window is taken as large as possible to 
increase statistics, its limit is due to the systematic error 
introduced because of the non-linear dependence of the 



counting rate on zenith angle (Alexandreas et al. 1993a). 
We performed a Monte Carlo (MC) calculation to find the 
largest radius that keeps the systematic error much lower 
than the statistical error due to the low counting rate. Re- 
sults are shown in Fig. |l]. The time intervals we chose are 
10 seconds, 1 minute and 4 minutes which cover the usual 
long duration GRBs. Dead time and low counting rate do 
not allow searches for GRBs with much shorter duration. 
The employed procedure has the advantage of being very 
sensitive but it is also quite time consuming. So in order 
to look for TeV-emission on longer time scales, a less sen- 
sitive search with a 1 hour time window and a classical 
non-overlapping rectangular grid of 0.5°x0.5° in celestial 
coordinates has been carried out. In this case the back- 
ground is estimated by mea ns of a MC method as de- 
scribed in Alexandreas et al. ( 1993a ). For every real event 
we generate 100 MC events with random directions follow- 
ing the acceptance function of AIROBICC (obtained with 
the events of a whole run, therefore a different run im- 
plies a different acceptance function) but with the same 
times as that of the original event. The number of MC 
events that fall within the ON window determines the 
background. 



Table 1. 90% confidence level upper limits for the inte- 
gral flux (in units of 10~ 8 cm~ 2 s _1 ) of any hypothetical 
GRB occurred in our data sample. They are tabulated 
depending on the time window and the declination band. 
Uncertainties are ^40%. 



Time scale 


10 sec 


1 min 4 min 


1 hour 


F UL (E >16TeV) 


4.5 


3.5 2.5 


0.35 


(9° < 5 < 49°) 








F UL (E >25TeV) 


2.8 


2.2 1.6 


0.32 


(-6° < 5 < 9°, 








49° < 5 < 64°) 









To obtain the significance for a hypothetical signal we use 
the probability distribution given in Alexandreas et al. 



(1993a) which is appropriate for low statistics (Poissonian 
regime). It is the probability of observing at least Non 
events in the source window, given the observed number 
of background events Noff, as a result of a background 
fluctuation: 



P(> Nqn I N OFF ) 
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noN 
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(new + NqffY- 
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where Non and Noff are the number of events in the 
source and background windows respectively. In case of 
random distribution of the events (i.e., no GRBs), a cu- 
mulative histogram of the number of trials with a chance 
probability lower than P as a function of -logio-P should 
follow a straight line with slope -1, which cuts the Y-axis 
at the height of the total number of windows (trials). A 
strong GRB or several weaker GRBs should therefore ap- 
pear as a deviation from the line. The significance of a de- 
viation can be calculated knowing that theoretically every 
bin content follows a binomial distribution. Actually this 
approach is only an approximation because the probability 
distribution is discrete and there is an oversampling (i.e., 
search windows overlap and thus trials are not indepen- 
dent) in the search with short time windows. Therefore we 
estimate the significance repeating the search over a MC 
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Table 2. List of GRBs observed in coincidence with satel- 
lites. 9 is the zenith angle of the GRB at the AIROBICC 
site. 



GRB 


Observer 


9 at trigger 


6 2 h later 


Coverage 


920525b 


BATSE 


10.0° 


31.4° 


Full 


920925c 


WATCH 


5.1° 


25.1° 


Full 


921118 


BATSE 


27.6° 


44.9° 


Partial 


930123 


BATSE 


29.6° 


10.0° 


Partial 



sample which is 10 times larger than the real one in the 
case of short time windows and 20 times larger in the case 
of 1 hour search window. 



4. Results 

The results for the all-sky search are shown in Fig. [| It 
presents the cumulative number of trials (search windows) 
against -logio-P for the all-sky search with the four time 
scales we have used. The results for the MC sample are 
also shown. Obviously, the experimental data set is con- 
sistent with MC expectations. The set of deviations has a 
chance probability greater than ~0.1. This result allows 
us to place an upper limit for the flux of any hypothetical 
GRB which may have occurred in the FOV during our 
observations. Depending on the declination band the re- 
sulting upper limit corresponds to two energy thresholds 
for every time scale. 

Assuming the same spectral index for the source and the 
CR flux and a steady emission during the time interval, we 
can estimate the int egral flux upper lim it with the formula 
(Karle et al. |l995b| ; Alexandreas et al. |l993b[) 



F<jl{E > E t h) 



Nul^on 



OFF 



1)0 



Fcr(E > E thhad ) 



where Fjjl is the 90% CL upper limit for the integral 
flux, E t h y and E t h had are the energy thresholds for 7 and 
hadrons respectively, Njjl is the 90% CL upper limit for 
the number of excess events in the source window as cal- 
culated by Aguilar-Bem'tez et al. ( 1992j ) and ft is the frac- 
tion of source events expected to fall in the source win- 
dow. Fcr is the known CR integral flux and is taken as 
1.8 • 10- 5 J B(TeV)~ 1 - 76±0 09 cm" 2 s" 1 sr" 1 (Alexandreas 



et al. 1993b; Burnett et al. 1990). We applied this calcu- 
lation to all ON-OFF windows and determined the high- 
est values. The resulting flux upper limits are listed in 
Table § Their uncertainties are estimated to be ~40% 
through comparison with a different method of flux calcu- 
lation. 

Additional constraints on the data sample as e.g. imposed 
by the burst detections of BATSE, WATCH and other 
GRB detectors in space, can serve as a tool to reduce the 
full data set and thus the expected statistical fluctuations. 
Hence, the sensitivity is enhanced allowing detection of 
weaker GRBs. We ther efore searched in the BATSE 3B 
catalog (Meegan et al. |l995[ ) and in t he WA TCH cata- 



Table 3. 90% confidence level upper limits for the inte- 
gral flux (in units of 10" 10 cm" 2 s" 1 ) of GRB 920525b in 
different time scales and for two scenarios: coincident and 
delayed emission. Uncertainties are ~40%. 



Time scale 


10 sec 1 min 


4 min 1 hour 


Ful{E >16TeV) 


89 8.9 


3.2 0.40 


(coincident, ±3 min) 






F VL (E >16TeV) 
(delayed, <2 h later) 


130 14 


3.2 







which were within the FOV of AIROBICC at the time 
they occurred or up to 2 hours later (because of possible 
delayed emission) . They are shown in Table 0. The sam- 
ple is then reduced to events within ±10° and ±3 minutes 
around GRB locations. We also looked for a delayed com- 
ponent for two hours after the initial triggers. The size of 
the search region in celestial coordinates has been chosen, 
a priori, large enough to have the uncertainties of GRB 
locations into account. The length in time of the search 
interval (for coincident and delayed emission) has been 
taken according to the observations in the GeV energy 



log (Castro-Tirado 1994; Sazonov et al. 1997) for triggers 



range (Hurley et al. 1994). 

The most significant excess in the full data set was found 
almost in coincidence with WATCH GRB 920925c, but 
from a direction 9° away from the most probable WATCH 
position. This observation is discussed in the next section. 
The results for the other three GRBs are shown in Fig. |[ 
It presents the cumulative number of trials against -logio-P 
for the four time scales used in the two search strategies 
(coincident and delayed emission) . The results for the MC 
sample are shown as a dashed line. No significant devia- 
tion is seen. This, again, allows us to place an upper limit 
for the integral flux of these GRBs. However, we only do 
so for GRB 920525b because it is the only one completely 
covered by the data set among the three GRBs. The upper 
limits for the integral flux are calculated at 90% CL in the 
same way as shown before. The results for the coincident 
and delayed emission with the four time scales used here 
are shown in Table ^. Their uncertainties are estimated to 
be ~40%. 

5. Observational results for WATCH GRB 
920925c 

As it was mentioned in the previous section, the search 
for emission from GRB 920925c had a surprising result. 
An excess was found searching for emission in coincidence 
with the WATCH trigger on all short time scales. No ex- 
cess has been found either in the search for delayed emis- 
sion or in the 1 hour time window. The results for the 
four time scales, in the two search modes (coincident and 
delayed emission), are shown in Fig. ^. Again the cumu- 
lative number of trials is plotted as a function of -logio-P 
and the dashed line shows the results for the MC sample. 
In this case significant deviations from the expectations 
are observed, especially in the 4 minute time scale (coin- 
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Fig. 3. Cumulative trials vs. probability in the search for the GRBs 920525b, 921118 and 930123. We show the results 
of the search for coincident and delayed emission. For GRB 930123 only the result for delayed emission is shown 
because AIROBICC started more than one hour after the burst trigger. The dashed line shows the MC results. No 
significant deviation appears. 
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Fig. 4. Cumulative number of trials (windows) versus probability in the search for GRB 920925c. We show the results 
of the search for coincident and delayed emission. The dashed line shows the MC results. A significant deviation 
appears in some of the plots. 



cident emission), which shows a deviation with a chance 
probability < 10 -4 . 

The search interval with the smallest chance probability 
(4 minutes time scale, coincident emission) is centered at 
a = 324.6° ± 0.3°, 5 = 16.8° ± 0.3° (J2000) and UT = 
22:45:21. Therefore it precedes the WATCH trigger by less 
than one minute and is about 9° away from the WATCH 
location. For this search interval we observe 11 events in 4 
minutes while 0.93 events are expected. The chance prob- 
ability for the background to yield such an excess, com- 
puted according to the formula shown above, is 7.7T0 -8 
and the Li and Ma significance (Li & Ma 1983) is 5.4<r 



(in both cases we take into account only 10 events inside 
the search interval, see above). This is the most signifi- 
cant excess seen in the whole data sample (compare Fig. |] 
with Fig. |2|). Correcting this probability with a trial fac- 
tor (using the MC) due to the search for the four GRBs 
in Table on several time scales and in a large solid angle 
region, yields a final probability of 3.3T0~ 3 (2.7a), i.e. the 
CL for this excess to be related to WATCH burst is 99.7% 
(neglecting the possibility of a second independent burst). 
The data registered with the HEGRA scintillator array 
does not show any excess, which may be due to its higher 
energy threshold and worse angular resolution. The time 
distribution of the events registered by AIROBICC in the 
position of the excess for the night of the GRB is plotted 
in Fig. |[ The distribution of the 11 events yielding the 
smallest chance probability is shown in detail in the inner 
part of the figure. Note that 7 out of the 11 burst events 
come within 22 seconds, 3 of which arrive within 0.25 sec- 



onds. 

Assuming that this excess is due to high energy emission 
we can calculate the integral flux in the same way as the 
upper limits, but replacing Njjl with Non~ cx{Noff + 1)- 
In this case the tentative mean integral flux above 16 TeV 
during the 4 minutes window is (9 ± 4)-10 -10 cm~ 2 s _1 . 
The burst as seen by the space detectors had a duration of 
~5 minutes exhibiting two main peaks. It was observed by 
WATCH and ULYSSES, thus reducing the possible loca- 
tions to an IPN (Interplanetary Network) annulus (Hurley 



1996). This annulus is obtained using the relative time of 



detection of both spacecrafts. The IPN ring passes through 
the WATCH 3er error circle and is 3° away from the excess 
observed with the AIROBICC array, see Fig. |^. 
Spectral data from WATCH in the interval 6-100 keV can 
be fitted by a power law spectrum. The fit yields a spec- 
tral index of 2.5±0.2 which, if naively extrapolated, does 
not predict any TeV flux detectable by AIROBICC. How- 
ever, the TeV spectrum may differ significantly from what 
one expects from the 100 keV extrapolation. Furthermore, 
the hypothetical emission at TeV energies precedes the 
WATCH observations and may therefore be due to an- 
other production mechanism. 

6. Discussion 

In order to add new clues that may clarify whether the ex- 
cess recorded with AIROBICC is related to GRB 920925c 



or not, the BATSE team (Kippcn 1996 ) investigated the 
reason why BATSE was not triggered by this GRB. There 
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Fig. 5. Counting rate of events in the position of the ex- 
cess detected by AIROBICC in coincidence with GRB 
920925c along the four hours of moonless night. A clear 
peak is seen at the time of the excess, the inner figure 
shows in detail the time distribution of the burst events. 
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Fig. 6. Map with the situation of the WATCH GRB 
920925c and the IPN ring as calculated from the WATCH 
and ULYSSES observations. The location of the excess 
detected by AIROBICC nearly coincident in time with 
WATCH is also shown. 



are two reasons why BATSE may have missed the burst: 
either it was not strong enough to overwrite the previous 
trigger (happened about one hour before) or the Earth (or 
the Moon) occulted the source. The BATSE team kindly 
provided us w ith a map of the sky, as seen by BATSE 
(Kippen |1996|) , at the time that GRB 920925c occurred. 





DEC \ \ \ 
30 \ \ 




........ \WATCH \ 

■ VoAIRObV 


:c 


RA 


30 





330 ': 


300 


1 


PNrec\ 










\— 


-30 j 


/ 













Fig. 7. Map with the situation of the WATCH GRB 
920925c and the AIROBICC excess. The recalculated IPN 
ring using the ULYSSES and AIROBICC times of detec- 
tion (see text) is also shown (IPNrec). 



It shows that the AIROBICC location was completely oc- 
culted by the Earth throughout the whole burst. Therefore 
the AIROBICC observation is compatible with BATSE 
not seeing the GRB. The WATCH position was on the 
Earth's limb at the beginning of the burst and was oc- 
culted about one minute later. The first peak of GRB 
920925c (which is within the first minute after the initial 
trigger) had lower peak flux than the previous BATSE 
trigger and thus it would have not produced a trigger over- 
write even neglecting the additional attenuation by the 
Earth's atmosphere. So the position given by WATCH is 
also consistent with the non-detection of the GRB with 
BATSE. 

A cross-check of burst positions given by BATSE, WATCH 
and IPN since 1991 to 1994 reveals disagreements for 
some events at a level of 5a or even more, so an error 
in the WATCH location cannot be discarded (if the er- 
ror is due to timing it would also affect the IPN loca- 
tion estimate). Furthermore it is interesting to note that 
assuming that the peak registered with AIROBICC cor- 
responds to the peak which triggered ULYSSES, we can 
recompute the IPN ring using the time of detection of 
AIROBICC (neglecting the WATCH- AIROBICC distance 
compared to the ULYSSES-Earth distance) yielding a new 
IPN ring which is compatible with AIROBICC but not 
with WATCH (see Fig. 0). 

At this moment it is not possible to decide the nature of 
the excess recorded with AIROBICC. Note that this GRB 
had a large fluence at keV energies, but was not very in- 
tense. The GRB occurred under a small zenith angle (12°) 
and could thus be studied with a low energy threshold of 
the AIROBICC array. The events of the excess differ sig- 
nificantly from background events as shown in Fig. 0. The 
figure plots the mean of the ratio of scintillator fired huts 
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Fig. 8. Difference between the events of the excess and 
background in the ratio of fired huts in the scintillator 
and AIROBICC arrays for groups of 11 events. The mean 
for the 11 events of the excess (vertical dashed line) is 
below the mean for background groups (Gaussian fit) at a 
level of 2.0(7. The implications of this result are discussed 
in the text. 

to AIROBICC fired huts for groups of 11 events. The his- 
togram fitted by a Gaussian function represents the groups 
of background events. The vertical dashed line, which is at 
2. Oct as given by the Gaussian fit, represents the mean for 
the 11 events of the excess. It has been shown (Arqueros 
et al. 1996) that the ratio particles to light for EAS at 
observation level is sensitive to the chemical composition 
and the difference observed in Fig. || points to a gamma 
origin of the excess. This is confirmed by a MC study in 
which simulated EAS are analyzed in the same way as 



real data (Martinez et al. 1995; Cortina 1997). The re- 
sults for real data and for MC are summarized in Table |]. 
If the excess has indeed gamma-ray nature its significance 
would be higher due to the much lower diffuse gamma- 
ray background compared to the charged particles back- 



ground (Karle et al. 1995c). On the other hand the angular 
separation between AIROBICC and WATCH locations is 
significant. 
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Table 4. Comparison between the ratio of scintillator 
fired huts to AIROBICC fired huts for groups of 11 events 
as predicted by the MC and as seen in the excess and in 
the background of the real data. The MC seems to confirm 
a gamma origin of the excess. 



<SCI/AIR> 


MC 




Data 






gammas 


hadrons 


excess 


background 


Mean 


1.26 


1.74 


1.29 


1.81 


RMS 


0.17 


0.23 




0.27 
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